Is my ODE a Painleve equation in disguise? 
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Abstract 

^ . Painleve equations belong to the class y"+a\ y' 3 +3a2 y' 2 +3a^ y'+a^ = 0, where a% = 

cii(x,y). This class of equations is invariant under the general point transformation 
x = &(X, Y), y = ^f(X, Y) and it is therefore very difficult to find out whether two 
equations in this class are related. We describe R. Liouville's theory of invariants 
that can be used to construct invariant characteristic expressions (syzygies) , and in 
' particular present such a characterization for Painleve equations I-IV. 

VO ' 

^ ■ 1 Introduction 

o 

Many phenomena in Nature are modeled by ordinary differential equations (ODEs). After 
such an equation is derived for some physical situation, the natural question is whether 
that ODE is well known, or at least transformable to a well known equation. For example, 
one would like to know if the equation is related to a known integrable equation, e.g, in the 
>- ! case of a second order ODE, to one of the Painleve equations, or to one of the equations 
in the Gambier|]/rrice[|2| list. 

Normally these lists contain only representative equations, e.g., up to some group of 
transformations. If the equation is of the form y" = f(x,y,y') the usually considered 
group of transformations is 

I = ^ (A) ' y = Ux)Y + uxy 

However, when we have an equation derived from some physical problem the required 
transformation may be more complicated. 

Here we consider the following class of equations 

y" + a\{x, y) y' 3 + 3a 2 (x, y) y' 2 + 3a 3 (x, y) y' + a 4 (x, y) = 0. (1) 

This class is invariant under a general point transformation 

x = $(X,Y), y = *(X,Y), (2) 
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we only need to assume that the transformation is nonsingular, i.e., 

A := <t> x ^ Y ~ ^y^x ^ 0. (3) 
The point transformation (0) prolongs to 

J = *x + *rY' (4) 

y" = [ [\& Y $ X - m x <& Y }Y" + [$ y fyy - ^y$yy]y' 3 

+ {§x^yy - 2$ X y^y + 2§ y ^xy - $yy^x)Y' 2 
+ (-^xx^y - 2^>xy^x + 2^ x ^xy + <$>y^xx)Y' 
+$x*xx-$xi%}/[^ + ^yf, (5) 

and when these are substituted into ([!]) the form of the equation (cubic in y') stays the 
same, but the coefficients ctj change as follows: 

Oi = ^($y^yy - $yyl£y + $ y A 4 + 3$ y ^yA 3 + 3$y^ y A 2 + ^yA{), (6) 
a 2 = ^[|($ X ^yy - 2$ X y^y + 2$y^ X y - $yy^ X ) + $ X $ y A 4 

+ (2$ x $y^y + ® Y ^x)A 3 + ($x^y + 2§ y ^x^y)A 2 + ^ X ^ Y A X ], (7) 

«3 = UK-^XxVy - 2$ X y^ X + 2$ X ^ X y + $y^ XX ) + $^$yA 4 

+ ($^y + 2$x®Y^x)A 3 + (2$ X ^ X ^y + <& Y ^ X )A 2 + ^^yA]], (8) 

a 4 = + + ® X A 4 + 3$^ X A 3 + m x ^ x A 2 + ^A). (9) 

Here Ai(X,Y) := aj($(X, K), Y)). Example: If we apply the transformation x = 
X + Y, y = XY to the equation y" = (in which a« = 0), then we find a[ = a 4 — 0, 
02 = 03 = i.e., the equation becomes K" + -^p(F' 2 + K') = 0. 

Since the coefficients of equation (TI]) transform in such a complicated way it is in 
general difficult to find a characterization of ([l|) that is invariant under a general point 
transformation (^). This classical problem was basically solved more than 100 years 
ago, the fundamental works being those by Liouville || in 1889, Tresse Q in 1894, with 
more modern formulations by Cartan || in 1924 and Thomsen M in 1930. Recent wave 
of interest on this classical problem started with 0. [For the restricted problem with 
x = $(X) see, e.g., §.] 



2 Relative invariants, absolute invariants and syzy- 
gies 

The invariants we are looking for must be constructed from the coefficients tij and their 
various derivatives. The transformation rules for aj were given in (|§-|9|) and we are now 
looking for some combinations that transform in a much simpler way. Let us consider 
some expression 

I[x, y] = I(a x , . . . , a 4 , d x a x , d x a ly d y at, d y a 4 , <9 2 ai, ...). 

Since = ai(x, y) this will be a function of x, y. Under the transformation (0) the ingre- 
dients ctj transform to 5i and the expression constructed from the transformed quantities 
in exactly the same way is 

I[X, Y] = 1(5]., . . . , o 4 , 9 X 5], . . . , <9 X 5], 5y5i, . . . , <9y5 4 , d x a x , ...). 
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This is now a function of X, Y and if it turns out that 

T[X,Y]=A n imx,Y),*(X,Y)], 

i.e., that I transforms, up to some overall factor, as by substitution then we say that 
J is a relative invariant of weight n. If furthermore n = or / = we say that / 
is an absolute invariant. Normally the weight of the relative invariant is indicated 
by a subscript. Later we will construct several sequences of relative invariants, e.g., 
^2n, n = 1,2,3,..., then using them we can construct a sequence of absolute invariants, 
in this case j 2n '■= i2n/i 2 - 

Although the absolute invariants transform by a simple substitution under (|2|) they are 
normally complicated functions of x, y and therefore in practice useless for classification. 
What we need are relationships between absolute invariants, i.e., syzygies. 

For example, it turns out later that for one particular equation — 6j4 + 4 = 0. 
This relationship is invariant under the transformation (0) and therefore it is an invariant 
characterization for that equation. Using this result we can say with certainty that any 
equation that does not satisfy this relationship cannot be transformed to first equation 
by any point transformation. The reverse is not true: several equations may satisfy the 
same syzygy. 

Note also that the syzygy polynomials can only have numerical coefficients, because 
parameter values appearing in the equation can be changed with the allowed transforma- 
tions. 



3 Construction of Invariants 

The construction of the relative and absolute invariants proceeds step by step as follows: 
First define 

n 22 = d^ — d y a 2 + 2(al — aia 3 ), (10) 
Ili2 = d x a 2 — d y a 3 + a 2 a 3 — aia 4 , (11) 
Tin = d x a 3 - d y a 4 + 2(aj - a 2 a 4 ) , (12) 



and then 



L 2 = d x Il 22 - d y Il 12 - a x Tl n + 2a 2 Ili2 - a 3 n 2 2, (13) 
L x = d x U 12 - d y Il n - a 2 U u + 2a 3 ni 2 - a 4 U 22 . (14) 

Now one finds that the Li transform according to 

%\ - a *A (LimX,Y)^(X,Y)]\ 

which is already a relatively simple transformation rule. From this we get the first result 
(R. Liouville): 

The property L\ = L 2 = is an absolute invariant and if it holds the equation can be 
transformed to Y" = 0. 
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If we define Oj and d x , d y to have weight |, then Lj are of weight |. Continuing with 
Lj and adding one derivative or Oj we find at weight | another pair 

Zi = -L Xy L x - 3L lx L 2 + 4:LiL 2x + 3a 2 L\ - Qa 3 L x L 2 + 3a 4 L2, (16) 
Z 2 = £2^2 + 3L 2y L 1 - 4L 2 L ly + 3a x L\ - Qa 2 L 1 L 2 + 3a 3 Ll, (17) 

that transforms similarly: 

Zl\ _ A 3 (*x (Z 1 MX,Y),^(X,Y)]\ 
%)' \®y VyJ \Z 2 [*(X,Y)MX,Y)]J 

Using these two we get the first semi-invariant 

v 5 = l[Z 1 L 2 -Z 2 L 1 ] 

= L 2 (L x d x L 2 - L 2 d x L x ) + L x {L 2 d y L x - L 1 d y L 2 ) 

-a x L\ + 3a 2 L\L 2 - 3a z L x h\ + a 4 L 3 2 , (19) 

which is of weight 5, i.e., transforms as u 5 = A 5 z/ 5 . 
Observation: For all Painleve equations v 5 — 0. 

Our ultimate aim is to provide an invariant characterization for the Gambier/Ince 
list of 50 equations. It contains some equations with L\ = L 2 = 0, and other equations 
besides Painleve's that have u 5 = 0, but it also contains many with i/ 5 ^ 0. [We note in 
passing, that in the standard form all Painleve equations also have the properties a± = 
and L 2 = 0, but these are not (semi) invariant characterizations.] 

Let us now continue with v 5 = 0, i.e. Z X L 2 = Z 2 L 1 . Let us define 

Ri = L x d x L 2 - L 2 d x L x + a 2 L\ - 2a 3 L x L 2 + a 4 L 2 2 , (20) 
w x = [LliUnU-UnL^ + RAiL^-Lld^ + L^asL.-a.L^/Li. (21) 

[If L x — there is a similar expression with L 2 as divider.] The expression w x is a 
semi-invariant of weight 1. 

Observation: For all Painleve equations w x — 0. 

The Gambier/Ince list contains also equations with v 5 = 0, w x ^ 0. 

We continue further with i/ 5 = 0, u>i = 0. A sequence of semi-invariants can now be 
constructed starting with 

i 2 = 2R 1 /L 1 + d x L 2 -d y L 1 , (22) 
with higher members given recursively by 

«2(m+i) = ^i dyi 2m - L 2 d x i 2m + 2mi 2m (d x L 2 - d y L x ). (23) 
If i 2 7^ a sequence of absolute invariants is given by 

32m = %2mi 2 m - (24) 

In the next section we will use the j 2m to characterize Painleve equations I-IV. 
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At this point we can see that the classification using j 2m cannot be sharp. Consider 
the special case y" + a^(x, y) = 0. Then Yl u = —d y a^ n 12 = IT22 = 0, L\ — dya i ,L 2 = 
0, z/ 5 = 0, Wi — 0, and 

i-2 = —dydi, «2(m+i) = L\ m+1 d y (i 2m / L 2 ™), 

The semi-invariants therefore depend only on d^a^ and its higher y-derivatives and are 
insensitive to the possible linear in y part in a^. 



4 Invariant characterization of Pj — Pjy 

The first steps were given above: All Painleve equations have the properties 1) at least 
one of Li, L 2 is nonzero, 2) z/ 5 = 0, 3) W\ = 0.Q If the candidate equation fails any of these 
properties it cannot be transformed to a Painleve equation using the transformation @. 
Here we want to go further and derive conditions which differentiate between Painleve 
equations. The following results have been derived using the symbolic algebra language 
REDUCE||. It should be noted that we give here only the lowest degree syzygy. 



4.1 Painleve I 

Any equation of the form 

y" = 6y 2 + fi(x)y + f (x), (25) 

has the property i 2 = 0. In principle i 2 is only a relative invariant, but since its value in 
this case is 0, this property is an absolute invariant. Equation (|25|) contains Painleve I 
for which fi = 0, fo = x. For all other Painleve equations i 2 ^ 0. 



4.2 Painleve II 

All equations of form 

y" = 2y 3 + f l (x)y + f (x), (26) 

have the property i 2 = 12, i± = 288 and therefore j'4 = 2 is the syzygy for this class of 
equations. [In fact it is easy to see that j2(m+i) = 2 m ml] Painleve II is contained as the 
special case fi = x, fo =const. 



4.3 Painleve III 

The following results are valid for Cases 12,13,13i in the Ince/Gambier classification, they 
all have four parameters, a, j3, 7, 5: 



/2 

y" - — --fy 3 -ay 2 -(3-5/y = 0, 

y 

y"- — + --jy 3 - -(m/ 2 + 0)- s/ y = 0, 

y x x 



y"- — - e x (ay 2 + (3)- e 2x ( 1V 3 4 */,,) = 0. 

y 



1 This was first observed by V. Dryuma in late 1980's 
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1: If any 3 of the parameters a, j3, 7, 8 are zero, then J2(m+i) — m ' 
2a: If 7 = 5 = or a = (3 = then 

j 6 -6j 4 + 4 = 0. (27) 

These two cases are connected: x = \X 2 , y = Y 2 takes (a, (3, 0, 0) — > (0, 0, a, /?). 

2b: If a = 7 = or (3 = 5 = then 

2j 8 - 13j 6 - 22 + 57j 4 - 21j 2 = 0. 

Transformation x = X,y = 1/Y takes (0, /3, 0, 5) -> 0, -5, 0). 

2c: If /3 = 7 = or a = 5 = then 

2j 8 -23j 6 -42 + 87j 4 -llj 2 = 0. 

Transformation x = X,y = 1/Y takes (a, 0, 0, 5) — > (0, —a, —5, 0). 

3: If a 2 <5 + /5 2 7 = (which contains 2a,2b) then we get a degree 12 relation 

4ji2 - 76jio - 1696 + 6840,74 - 2600 j 6 - 5640 j\ 

+590j 8 + 2220j 4j6 + 450j 4 3 - 155j 4 j 8 - 115j 2 = 0. 

4: For the generic case we obtain 

4j 14 - H2ji2 - 231j 10 j 4 + 1274j 10 - 385j 8 j 6 + 4795j 8 j 4 - 7910j 8 + 3255j 2 

+3570j 6 j^ - 39060j 6 j 4 + 30240j 6 - 15330j 4 3 + 78120j^ - 71736j 4 + 15264 = 0. 

4.4 Painleve IV 

Let us again consider a more general form 

13 1 
y" = ^y' 2 + e ^y 3 + 4 ( e 2^ + e 3)v 2 + 2(e 4 x 2 + e 5 x + e 6 )y + e 7 -, 
2y 2 y 

Painleve IV corresponds to e\ = e 2 = e 4 = 1, e 3 = e$ = 0, = —a, e 7 = —/3 2 /2. The 
classification is sensitive only to e±, e 2 , e 3 , e 7 . 

1: If e 1 = e 2 = e 3 = e 7 = then L\ = L 2 = 

2a: If e 2 = e 3 = e 7 = then j 2 (m+i) = (I)™ 7 ™' 

2b: If e 2 = e 3 = ei = then j 2 ( m +i) = {fT m ^ 

2c: If ei = e 7 = then j 2 (m+i) = 2 m m! 

3a: If e 7 = then 

j 8 -llj 6 - 24 + 46j 4 -7j 2 = 0. 
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3b: If ei = then 



5j 8 - 95j 6 - 136 + 286j 4 + 21jj = 0. 

3c: if ti = e% = then 

8125j 10 - 2165j 8 - 2495j 6 j 4 + 13096j 6 + 12872 j\ - 36160 j 4 + 11904 = 0. 

4: Generic case: 

3ji4 + 1955j 10 - 277jioj 4 - 120ji 2 + 31488 - 131008j 4 + 57952j 6 + 103600j'| 

-15062j 8 - 40486j 4 j6 - 18772j 4 3 + 5236j 4 J 8 + 990j 6 2 + 3775j 4 2 j 6 - 75j 6 j 8 = 0. 

5 Conclusions 

We have presented here the first results of a project aiming to derive syzygies for every 
equation in the Gambier/Ince list. The biggest problem in finding the characteristic 
expressions is that we have to solve huge sets of rather complicated nonlinear algebraic 
equations. For Painleve V the generic expression is of degree higher than 26 (which 
is as far as we checked) and for Painleve VI presumably still much higher. When this 
classification is eventually finished: 

• We get an algorithmic method for finding out if a given equation of type (|l|) has 
any change of being transformed into one of the equations in the Gambier/Ince list. 

• We also get a classification of the equations into essentially different sub-cases (c.f. 

The first results that have already been obtained suggest some interesting open ques- 
tions: 

• Is the type of "integrability" different in the classes a) v 5 = wi = 0, b) u 5 = 0, W\ ^ 
0, and c) z/ 5 ^ 0? 

• What does the degree of the minimal syzygy tell us about integrability? 
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